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COMBINATIONS  CF  TEMPERATURE  AND 
AXIAL  COMPRESSION  REQUIRED  FOR  BUCKLING 
OF  A  RING-STIFFENED  CYLINDER 
By  Melvin  S.  Anderson 


SUMMARY 


A  theory  is  presented  to  predict  the  buckling  temperature  of  an 
axially  compressed,  uniformly  heated  ring-stiffened  cylinder.  The 
cylinder  buckles  because  of  the  interaction  of  the  axial  stress  due  to 
applied  compressive  loads  and  the  circumferential  stress  resulting 
from  restraint  of  thermal  expansion  by  the  rings.  Buckling  charts 
covering  a  vide  range  of  cylinder  proportions  are  presented  for  both 
clamped  and  simply  supported  cylinders.  The  buckling  temperature  for  a 
given  axial  loading  is  determined  from  a  simple  equation  involving  a 
coefficient  given  in  the  buckling  charts  and  the  radius-thickness  ratio 
of  the  cylinder. 


INTRODUCTION 


One  aspect  of  the  structural -aerodynamic  heating  problem  is  the 
effect  of  thermal  stresses  on  buckling.  Reference  1  reports  the  results 
of  an  experimental  investigation  which  indicated  the  effect  of  thermal 
stresses  on  buckling  and  maximum  strength  for  a  variety  of  structures. 
The  results  for  the  ring-stiffened  cylinders  indicated  the  compressive 
buckling  stress  was  reduced  by  thermal  stresses.  These  thermal  stresses 
were  in  the  circumferential  direction  and  caused  by  the  restraint  of 
thermal  expansion  in  the  vicinity  of  the  rings. 

Reference  2  presented  a  theoretical  analysis  of  the  buckling  of  a 
simply  supported  cylinder  due  to  this  type  of  thermal  stress  alone.  The 
calculated  buckling  temperature  of  a  steel  cylinder  with  a  radius - 
thickness  ratio  of  about  5OQ  "was  found  to  be  over  2,500°  F)  thus,  there 
appeared  to  be  no  problem  in  the  practical  temperature  range  of  the 
material.  Calculations  for  other  proportions  were  not  given.  For 
larger  values  of  radius -thickness  ratio,  the  buckling  temperature  is 
lower  and,  as  shown  later,  thermal  buckling  can  occur  for  cylinders  of 
practical  proportions.  The  present  paper  presents  a  theoretical  analysis 


for  buckling  of  a  cylinder  subject  to  axial  compression  as  well  as  thermal 
stress  for  both  simply  supported  and  clamped  end  conditions.  A  comparison 
of  the  theory  with  the  exi)erlmental  results  of  reference  1  is  also  given. 
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cross-sectional  area  of  ring 
coefficient  in  deflection  function 
constant  (see  eq.  (6)) 

plate  flexural  stiffness, 

modulus  of  elasticity 

circumferential  stress  coefficient  for  arbitrary  stress 
distribution 


temperature  coefficient. 


12  (oT  -  ai^TR)(l  -  d2)/L'.2 


^  ) 


modified  temperature  coefficient  applicable  to  cylinders  'with 
flexible  rings 


axial  stress  coefficient. 


12ax  1  -  |i‘ 


clrc'umferential  stress  coefficient  for  uniform  stress  distri- 

^y( 


bution. 


jt^E 


length  of  cylinder  between  rings 


=  [(m  -  1)2  + 


+  (a  -  1)^ 

[(m  -  1)2  +  2 


(m  -  l)2k2j 


m,n,l 


Integers 
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Q  operator  defined  in  equation  (A5) 

r  radius  of  cylinder 

8jj  Fourier  comi)onents  of  circumferential  stress 

T  temperature 

temperature  difference  between  cylinder  wall  and  ring,  T  -  Tj^ 

t  thickness  of  cylinder  wall 

Vjjj  deflection  ftinction 

w  radial  deflection  of  cylinder  wall 

x,y  axial  and  circumferential  coordinates  (see  fig.  12) 

Z  cyli  nder  curvature  parameter,  ~  v/l  - 

rt  ’ 

a  coefficient  of  thermal  expansion 

p  ratio  of  cylinder  length  to  circumferential  buckle  length,  L/X 

7  axial  stress  ratio)  ratio  of  axled,  stress  to  the  classical 

Or 

cylinder  buckling  stress,  - - - 

^5(1  -  n^h 

7'  modified  axial  stress  ratio  based  on  experimentally  observed 

reduction  in  cylinder  buckling  stress  from  classical  value 

6jnn  Kroenecker  delta,  =  0  \dien  m  jii  n  and  =  1  when 

m  =  n 

6  =  fifz  (1^ 

0  =  Jl  +  7 

circumferential  buckle  length 

H  Poisson's  ratio  (p  eqml  O.3  was  used  for  all  nvonerical 

calcxilatlons) 
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ffjj,  axial  coopressive  stress 

Oy  circumferential  compressive  stress 

dy  average  circumferential  compressive  stress 


T 

T' 

i 


temperat\ire  coefficient^ 


(»!  -  aElB)E 


modified  temperature  coefficient  applicable  to  cylinders  vitii 
flexible  rings 

deflectional  stiffness  of  ring  in  radial  direction 


differential  operator, 

dx^  dx^dy^ 

inverse  operator  defined  such  that  'v^\Aw  =  V 
Subscripts*, 
max  maximum 

R  refers  to  ring 

0  evaluated  at  t  eq,ual  zero 

Bars  over  symbols  indicate  average  values. 
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THEORY  AND  ASSUMPTIONS 


The  cylinder  is  assumed  to  have  many  rings  vhich  divide  the  cylinder 
into  bays  of  length  L.  The  cylinder  vail  is  at  a  uniform  temperatxire 
and  the  rings  are  at  some  lover  temperature.  This  difference  in  tempera¬ 
ture  causes  circumferential  stresses  vhich  vary  in  the  axial  direction 
and  have  the  largest  magnitudes  in  the  vicinity  of  the  rings.  Combina¬ 
tions  of  axial  compression  and  temperature  (circumferential  thermal 
stress)  necesseury  to  cause  buckling  vere  determined  from  the  analyses 
presented  in  the  appendixes. 

In  appendix  A,  the  problem  considered  is  buckling  of  a  cylinder 
subjected  to  a  uniform  axial  stress  and  a  circumferential  stress  uhat 
varies  in  the  axial  direction  and  is  expressed  as  a  general  Fourier 
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series.  In  appendix  B,  the  circumferential  thermal  stress  resulting 
from  the  ten^jerature  difference  between  the  cylinder  wall  and  ririga  is 
determined}  this  result  is  expressed  as  a  Fourier  series  in  order  to  be 
applicable  to  the  buckling  analysis  of  appendix  A.  One  of  the  major 
assumptions  of  the  whole  analysis  is  in  applying  the  thermal  stress 
results  to  the  buckling  problem.  It  is  assumed  that  circumferential 
stresses  produced  by  the  teii5>erature  rise  of  the  cylinder  can  be  used 
in  the  buckling  analysis.  However,  the  deflections  caused  by  the  tem¬ 
perature  rise  are  neglected  Inasmuch  as  a  perfect  cylinder  is  assumed 
in  the  buckling  analysis.  Some  support  of  this  assumption  is  given  by 
the  fact  that  these  deflections  ore  generally  small  and  not  compatible 
with  the  final  buckle  pattern. 

In  the  stability  analysis  of  appendix  A,  the  modified  equation  of 
equilibrium  presented  by  Batdorf  (ref.  5)  was  used.  The  buckle  deflec¬ 
tion  was  expressed  as  a  trigonometric  series  and  an  infinite  stability 
determinant  resulted  when  the  series  was  substituted  into  the  equilibrium 
equLation.  For  simply  supported  cylinders,  the  results  obtained  using  the 
modified  equation  are  the  same  as  those  that  would  be  obtained  by  using 
Donnell's  eighth-order  differential  equation,  However,  as  pointed  out  in 
reference  5,  the  use  of  Donnell's  equation  for  clamped  cylinders  may  lead 
to  divergent  trigonometric  series  but  this  problem  is  not  encountered 
with  the  use  of  the  modified  equation. 


The  usxial  assumptions  of  small-deflection  cylinder  theory  together 
with  the  assumptions  implied  by  the  use  of  the  modified  equilibrium 
equation  are  present  in  the  analysis.  These  assumptions  include  cer¬ 
tain  restrictions  on  the  in-plane  displacements  at  the  end  of  the 
cylinder.  For  the  deflection  functioriR  used  in  this  paper,  the  following 
conditions  are  implied  on  the  in-plane  displacements  at  the  ends  of  the 
cylinder;  for  simply  supported  cylinders,  free  displacement  in  the  axial 
direction  and  zero  displacement  in  the  circumferential  direction;  and  for 
clamped  cylinders,  zero  displacement  in  the  axial  direction  and  free  dis¬ 
placement  in  the  circumferential  direction.  The  condition  of  free  dls- 


I  o  r 


not  seem  to  be  too  realistic  but  it  is  pointed  out  in  reference  5  that 
there  is  little  difference  in  critical  stress  obtained  with  this 
boundary  condition  and  the  critical  stress  obtained  by  other  investi¬ 
gators  using  a  boundary  condition  that  all  edge  displacements  were  zero. 


The  rings  are  assumed  to  be  rigid  against  radial  loads  but  are 
allowed  to  expand  because  of  a  temperature  rise .  A  correction  which 
allows  for  the  ring  expansion  that  occurs  when  the  cylinder  is  heated  and 
loaded  is  given  in  a  later  section. 
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THEORETICAL  RESULTS 


The  results  of  the  analysis  are  presented  in  figures  1  and  2  and  in 
tab\ilar  form  in  table  I .  Results  are  presented  as  a  buckling  temperature 
coefficient  t  plotted  against  the  cylinder  curvature  parameter  Z  for 
various  values  of  y ,  the  ratio  of  applied  axial  stress  to  the  so-called 
classical  buckling  stress.  The  buckling  temperature  can  '^e  calculated 
as 


T 


=  5!Et„  +  1 1 

a  ^  a  r 


or^  if  the  coefficient  of  thermal  expansion  of  ring  and  skin  can  be  con¬ 


sidered  equal 


then 


AT  =  i  ^ 
a  r 


(2) 


In  practice,  T  vill  represent  the  average  skin  temperature,  Tj^  the 
average  ring  temperature,  and  AT  the  difference  between  the.se  average 
temperatures.  The  values  of  t  shown  are  the  lover  of  the  two  values 
foxind  for  the  symmetrical  and  anti symmetrical  deflection  patterns  and 
have  been  minimized  with  respect  to  the  wavelength  parameter  p,  the 
ratio  of  the  cylinder  length  to  circvunferential  buckle  length.  Although 
the  precise  minimum  was  not  found,  the  increments  in  P  were  suitably 
small  so  that  little  error  is  present  in  using  the  lowest  value  that  was 
calculated. 

Inasmuch  as  the  thermal  stress  is  concentrated  over  the  rings,  many 
terms  in  the  deflection  function  are  required  to  describe  the  buckle 
pattern  accurately,  especially  at  large  values  of  Z.  Calculations  were 
made  on  an  IBM  JO9O  electronic  data  processing  system  and  by  increasing 
the  order  of  the  stability  determinant  suitable  accuracy  was  achieved 
except  in  some  cases  as  indicated  later.  The  results  given  are  such 
that  using  a  determinant  of  one  less  order  would  yield  a  value  of  t  no 
more  than  5  percent  greater .  However,  except  for  the  larger  values  of 
Z,  most  of  the  values  of  t  given  in  the  table  differ  by  less  than 
1  percent  from  values  obtained  by  using  a  determinant  of  one  less  order. 
The  order  of  determinant  used  as  well  as  the  value  of  P  is  given  in 
the  table .  In  some  cases  it  was  not  possible  to  get  an  accurate  value 
of  T  as  the  size  of  the  determinant  was  limited  by  the  computer  program. 
These  points  are  omitted  from  the  table  but  extrapolations  based  on 
established  calculations  are  shown  by  dashed  lines  in  the  figures.  At 
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low  values  of  Z  where  huckling  is  associated  with  plate  action  rather 
than  cylinder  action,  it  is  possible  for  7  to  be  greater  than  1  and 
calciilations  for  these  points  are  given  in  the  table  only.  The  results 
given  in  figures  1  and  2  show  that,  as  Z  is  increased,  7  being  held 
constant,  t  becomes  independent  of  Z.  This  is  true  for  all  Values  of 
7  for  the  simply  supported  cylinder.  However,  it  is  not  completely 
definite  that  this  is  true  for  the  smaller  values  of  7  in  the  case  of 
the  clamped  cylinder  because  of  the  difficulty  in  obtairjing  an  accurate 
value  of  T.  Although  values  are  given  only  to  Z  =  1,000,  calculations 
for  Z  greater  than  1,000  indicate  the  trend  shown  would  continue. 

The  value  of  t  for  7  equal  to  zero  can  be  used  to  determine  the 
temperature  at  which  a  cylinder  will  buckle  in  the  absence  of  aocial 
stress.  For  the  range  of  Z  in  which  t  is  independent  of  Z,  the 
buckling  temperat-ure  is  a  function  of  r/t  and  the  coefficient  of  thermal 
expansion  a,  as  indicated  in  figure  5  where  the  temperature  rise  neces¬ 
sary  to  cause  buckling  of  aluminum  and  stainless-steex  cylinders  is 
plotted  against  r/t.  A  constant  value  of  a  was  assumed  for  ring  and 
skin  in  each  case.  At  low  values  of  r/t  the  buckling  temperature  is 
very  high  and  is  beyond  the  useful  temperature  range  for  the  material. 
However,  at  high  values  of  r/t  the  buckling  temperature  can  be  low 
enough  that  thermal  buckling  would  have  to  be  considered  in  the  design  of 
a  heated  cylinder.  It  sho\ild  be  noted  that  the  curves  in  figure  5  also 
apply  to  values  of  7  up  to  0.7. 

The  buckling  temperature  of  a  simply  supported  cylinder  subjected 
to  heating  alone  can  also  be  obtained  from  the  analysis  of  reference  2. 

In  the  numerical  example  of  reference  2,  a  buckling  temperature  rise  of 
2,580°  F  was  calculated  for  a  steel  cylinder  with  r/t  equal  to  5OO. 

It  can  be  seen  that  this  temperature  is  about  three  times  the  value  shown 
in  figijire  5*  The  difference  is  probably  due  to  the  boiondary  conditions 
used  in  calc\ilating  the  thermal  stress  as  the  other  parts  of  the  two 
analyses  are  identical.  Reference  2  assumed  a  one-bay  cylinder  and  cal¬ 
culated  the  thermal  stress  for  simply  supported  ends.  The  present 
analysis  assumes  many  bays;  this  assvimption  effectively  clamps  the  ends 
as  far  as  thermal  stress  is  concerned.  (See  appendix  B.) 

For  the  range  of  Z  in  which  t  is  independent  of  Z,  the  results 
given  in  figures  1  and  2  can  be  cross  plotted  to  give  curves  for  the 
variation  of  t  with  7  as  shown  in  figure  4.  As  expected,  when  axial 
compression  is  the  dominant  factor  in  buckling,  there  is  little  difference 
between  a  clamped  and  simply  supported  cylinder.  However,  at  small  values 
of  7,  circumferential  stress  is  the  most  important  factor  in  buckling  and 
there  is  a  significant  difference  between  the  results  for  a  clamped 
cylinder  and  the  resxilts  for  a  siit^jly  supported  cylinder.  This  difference 
is  consistent  with  that  obtained  for  xmifonn  circ\miferential  compression. 
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The  curves  of  figure  4  indicate  that  an  appreciable  amount  of  com-  ^ 

pression  can  be  applied  without  reducing  the  ten^jerature  at  which  the 
cylinder  will  buckle.  For  the  clamped  cylinder,  the  buckling  temperature 
can  actually  Increase  with  applied  load-  This  result  appears  to  be  con¬ 
trary  to  intuition  but  it  can  be  shown  that  the  average  circumferential 
stress  is  declining  even  though  the  temperature  is  increasing. 

It  is  of  interest  to  note  that  in  figure  4  t  is  zfero  for  7 
approximately  eq.ual  to  0.95*  This  result  means  that  the  theoretical 
buckling  stress  for  an  unheated  cylinder  (t  =  O)  is  about  95  percent  of 
the  classical  value  because  of  the  restraint  of  Poisson's  expansion  at 
the  rings  or  bulkheads. 

The  sharp  break  in  the  curves  of  figure  4  is  the  result  of  a  mode 
change  which  is  illustrated  in  figure  5*  The  coefficient  t  is  plotted 
against  3  for  a  clampied  cylinder  with  Z  equal  to  50  and  7  equal  to 
0.6  and  0.7.  For  each  value  of  7,  curves  corresponding  to  the  two  ♦ 

lowest  buckling  modes  are  shown.  The  controlling  mode  for  7  equal  to 
0.6  has  a  minimim  value  of  t  at  p  equal  to  5*4.  For  7  equal  O.7 

there  is  another  minimum  in  the  same  region;  however,  there  is  a  lower 

minimum  at  p  equal  to  2.4.  It  can  be  seen  that  a  mode  change  has  taken 
place  and  that  a  sharp  break  will  occur  in  a  plot  of  the  minim-urn  value 
of  T  against  7  when  going  from  one  mode  to  the  other.  At  the  lower 
values  of  7,  where  buckling  is  due  primarily  to  circumferential  stress, 
the  T  -  p  curves  for  the  controlling  mode  shape  are  characterized  by 
rather  flat  minimuns  and  require  large  stability  determinants  for  con¬ 
vergence.  Calc-ulations  of  the  associated  buckle  deflection  indicate 
large  deflections  near  the  ends  where  the  thermal  stress  is  concentrated 
and  essentially  zero  deflection  a-way  from  the  ends.  At  larger  values  of 
7,  the  mode  shape  changes,  and  the  t  -  p  curves  have  sharper  minimums 
and  do  not  require  as  large  a  stability  determinant  as  the  lower  values 
of  7.  For  this  case  the  calculated  buckle  deflections  are  in  a  fairly 

regilai-  -wave  pattern  over  the  entire  length  of  the  cylinder  as  is  the 

case  for  compressive  loads  acting  alone. 


Effect  of  Thermal  Stress  Distribution  on  Buckling 

As  mentioned  previously,  the  thermal  stress  is  concentrated  in  the 
vicinity  of  the  rings  as  indicated  in  figure  6.  The  circumferential 
stress  at  a  ring  (x  =  0)  can  be  obtained  from  equations  (AJ)  and  (B5)  as 


byCo)  =  +  oET  -  opETp 


(5) 


tH  O^VD  O 
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circimferential  stress  distritution  was  ottained  from  appendix  B 
and  is  shown  for  two  values  of  Z  and  two  values  of  y  in  figure  6  as 
the  ratio  of  circumferential  stress  to  circumferentisil  stress  at  a  ring 
plotted  against  x/L.  A  study  of  figure  6  and  additional  thermal  stress 
calcvilatlons  indicate  that  for  cylinders  of  a  given  radius  and  thickness 
the  thermal  stress  distribution  in  the  vicinity  of  the  rings  is  xanaffected 
by  increases  in  length  beyond  a  certain  value,  and  the  thermal  stress 
away  from  the  rings  is  essentially  zero  for  these  "long"  cylinders.  This 
behavior  is  probably  the  reason  that  the  buckling  temperatxire  is  independ¬ 
ent  of  length  for  the  larger  values  of  Z. 


The  average  value  of  the  circumferential  stress  Oy  is  shown  in 


figure  6  as  the  line  denoted 


It  is  of  interest  to  compare  the 


average  value  of  the  thermal  buckling  stress  with  the  uniform  circum¬ 
ferential  stress  that  causes  buckling.  Values  of  ky,  the  buckling 


coefficient  for  uniform  circumferential  compression,  can  be  obtained  as 
a  special  case  from  the  analysis  in  appendix  A.  Inasmuch  as  the  thermal 
stress  over  the  ring  can  be  obtained  from  the  buckling  coefficient  k/p. 


a  buckling  coefficient  corresponding  to  the  average  thermal  stress  is 


A  comparison  of  the  buckling  coefficient  for  uniform  circum' 


ferential  compression  ky  with  the  buckling  coefficient  for  the  average 
thermal  stress  — 7— is  given  in  figure  7>  ^  buckling  coefficient 

Oy(0) 

based  on  cylinder  length  was  used)  this  manner  of  presentation  corre¬ 
sponds  to  the  way  in  which  results  are  usually  presented  for  uniform 
circumferential  compression. 


At  low  values  of  Z  the  thermal  stress  distribution  is  almost 
uniform  so  ^ ^ j k^i  and  ky  are  armost  idenhicai  uuu.  approacii  ciic  vu.j_u^ 

y 

for  a  flat  plate  as  Z  approaches  zero.  As  Z  increases,  the  thermal 
stress  distribution  becomes  similar  to  that  shown  in  figure  6  and  it  can 
be  seen  from  figure  7  that  the  average  thermal  stress  at  buckling  is 
about  three  to  four  times  the  buckling  stress  for  a  uniform  distribution 
for  both  simply  supported  and  clamped  cylinders. 


Effect  of  Ring  Flexibility  on  Buckling 

The  results  given  in  table  I  and  figures  1  and  2  apply  to  rings 
which  do  not  deflect  radially  when  the  skin  is  heated  or  loaded  in  axial 
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compression.  The  deflection  ein  actual  ring  undergoes  can  he  taken  into 
account  as  indicated  in  appendix  B.  The  buckling  temperature  coeffi¬ 
cient  t'  for  nonrigid  rings  is  given  in  terms  of  the  rigid-ring  solu¬ 
tion  plus  an  additional  term  involving  the  average  circumferential 
stress  in  the  skin  and  the  ratio  of  skin  and  ring  stiffnesses.  The 
coefficient  t'  is  given  as 


=  T  + 


^3(1  -  ErAr 


ELt 


(4) 


It  may  he  necessary  in  applying  eq,uatlon  (^4-)  to  use  an  effective  ring 
area  Aj^  if  the  entire  ring  does  not  fully  participate  in  resisting  the 


expansion  of  the  skin.  Values  of 


By(0) 


which  are  necessary  for  the 


use  of  equation  (4)  have  been  calculated  and  are  shown  in  figure  8 

5v 

plotted  against  Z,  It  is  seen  that  for  Z  smallj,  —  approaches  1 
because  the  stress  distribution  is  nearly  uniform.  As  Z  increases, 

^  becomes  proportional  to  l/y^  which  can  be  considered  as  propor¬ 
tional  to  1/L  for  r/t  constant.  Such  a  result  conforms  to  the 
observations  made  in  the  previous  section,  that  the  thermal  stresses  are 
concentrated  near  the  rings  and  are  independent  of  cylinder  length  for 
the  larger  values  of  Z.  Thus  for  a  given  radius  and  thickness,  the  net 
compression  force  resulting  from  the  thermal  stress  is  constant  and  the 
average  value  of  thermal  stress  is  proportional  to  1/l  as  indicated  in 

n  ^V 

figure  8.  The  oscillations  of  — at  higher  values  of  7  are  about 


Oy(C» 


I  i.ii  i.  ^ 


one-half  cycle  out  ox  phase  with  the  uBcxiiatiouD  oi 
and  2.  Thus,  t — which  may  be  thought  of  as  an  average  stress 

Oy(0) 

buckling  coefficient,  is  a  smoother  function  of  Z  than  t.  In  the 
region  where  — is  proportional  to  1 //Z,  equation  (4)  can  be 
written  Independent  of  length  for  intermediate  to  large  values  of  Z  as 


T' 


T  + 


ayi  -  7  £t5  _E. 

V  %  % 


(5) 
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APPLICATION  OP  RESUI/PB  AND  COMPARISON  WITH  EXPERIMENT 


Modified  Theory 

Because  of  the  discrepancy  hetween  experimental  results  and  small- 
deflection  buckling  theory  for  cylinders  under  axial  compression  the 
present  theory  cannot  be  expected  to  predict  experimental  results 
accurately  for  cases  in  vhich  axial  compression  is  the  predominant  factor 
in  buckling.  It  is  necessary  therefore  to  modify  the  res\alts  similar  to 
what  has  been  done  for  cylinders  under  compression.  Of  course,  any 
modification  should  be  such  that  experimental  data  for  cylinders  -under 
uniform  compression  agree  with  the  modified  theory  -when  t  equals  zero. 
Reference  presents  an  empirical  approach  to  obtain  a  reduced  buckling 
stress  for  cylinders  under  uniform  compression.  The  buckling  stress  is 
given  by 


In  equation  (6) ,  B  is  an  experimentally  determined  constant  which  may 
be  thought  of  as  a  measure  of  initial  imperfections  present  in  the 
cylinder  wallj  if  B  is  zero,  equation  (6)  is  the  usual  small-deflection 
buckling  equation.  If  a  prime  is  used  to  indicate  a  modified  buckling 
ratio  and  a  subscript  zero  used  to  indicate  t  being  equal  to  zero, 
equation  (6)  can  be  written 


7 


1 

1  +  6b^ 


(7) 


A  simple  extension  of  this  eqxiatlon  that  will  give  a  reduced  value  of  y 
for  any  value  of  t  is 


7' 


(8) 


From  the  shape  of  the  t  -  7  interaction  c\irve  given  in  figure  4  it  can 
be  seen  that  equation  (8)  wi3J.  pro-vlde  the  proper  amount  of  correction  in 


12 


# 


the  vicinity  of  t  equal  zero  hut,  at  lower  values  of  y  where  the 
curve  is  relatively  flat  and  the  circxnnferential  stress  is  the  pre¬ 
dominant  factor  in  buckling,  little  change  in  the  interaction  curve 
will  he  obtained.  Inasmuch  as  the  buckling  stress  of  a  cylinder  loaded 
in  circumferential  con^iressicu  agrees  fairly  well  with  theory,  a  correc¬ 
tion  should  not  be  required  where  circumferential  compression  is  the 
major  factor  in  buckling.  Therefore,  equation  (8)  may  be  expected  to 
give  a  reasonable  estimate  of  the  compressive  buckling  coefficient  when 
a  cylinder  is  heated  despite  the  limited  amount  of  experimental  data 
presently  available  to  substantiate  such  a  correction.  However,  seme 
Judgment  must  be  exercised  in  applying  equation  (8)  dependent  on  the 
shape  of  the  interaction  curve.  For  example,  in  a  region  of  the  curve 
where  x  Increases  with  y,  as  occurs  for  clamped  cylinders,  the  type 
of  correction  required  is  not  clear. 


Comparison  of  Theory  and  Experiment 

In  reference  1,  buckling  test  results  are  given  for  several 
2024— T3  aluminum  cylinders  loaded  in  bending  and  uniformly  heated. 

Although  the  theory  of  appendix  A  applies  to  uniform  compression  only, 
it  is  assumed  that  bending  tests  can  be  compared  directly  with  the  theory 
for  uniform  compression.  Theoretical  studies  such  as  reference  5  have 
indicated  that  the  maximum  compressive  stress  at  buckling  for  a  cylinder 
loaded  in  pure  bending  is  essentially  the  classical  buckling  stress  for 
uniform  compression.  Experimental  results  have  shown  a  small  difference 
between  compression  and  bending  (see  ref,  5)^  ^ut  this  is  accounted  for 
by  the  choice  of  B  in  equation  (8) .  A  comparison  of  the  results  of  ref¬ 
erence  1  with  the  theory  of  the  present  paper,  including  the  reduction  in 
axial  buckling  stress  Indicated  in  equation  (8),  is  shown  in  figure 
The  extreme  fiber  corapreesive  stress  at  buckling  is  plotted  against  the 
maximum  cylinder  temperature.  The  theoretical  curves  were  calculated  by 
assuming  singly  supported  ends,  although  assuming  clamped  ends  would 
yield  essentially  the  same  results  for  the  range  of  variables  encountered. 
The  correction  for  ring  ij.exibii.lty  ^eq.  '\^j}  was  ut>cu  out  axu  noi.  euLUSi' 
the  results  appreciably  from  the  rigid-ring  case. 

In  order  to  calculate  the  curves  shown  in  figure  9/  it  is  necessary 
to  know  the  variation  of  the  physical  quantities  E  and  a  with  temper¬ 
ature  as  well  as  the  relationship  of  the  skin  and  ring  average  tempera¬ 
tures  to  the  maximum  cylinder  temperature.  The  values  of  a  and  E 
that  were  used  in  the  calculations  are  shown  in  figure  10.  The  difference 
in  expansion  coefficient  for  the  skin  and  the  lower  temperature  ring  was 
neglected)  therefore,  only  the  difference  in  average  ring  and  skin  temper- 
attire  is  required  which  is  given  in  figure  11  as  a  function  of  maximum 
cylinder  temperature.  It  was  not  necesseiry  to  use  equation  (8)  to  modify 
the  theory  for  the  cylinder  having  L/r  equal  to  l/4  inasmuch  as  little 
difference  between  small-deflection  theory  and  experiment  was  found  for 
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this  proportion  as  indicated  by  the  room- temperature  test  point.  For 
the  remaining  cylinders  a  value  of  B  equal  O.OOOJO  was  used  in  equa¬ 
tion  (8) .  It  was  found  in  reference  ^  that  this  value  gives  results 
which  are  in  the  middle  of  the  experimental  scatter  band  for 
aluminum  cylinders.  The  point  for  L/r  equal  to  1  which  is  below  the 
theory  would  lie  above  the  theory  if  the  lower  limit  of,  the  scatter  band 
of  the  data  of  reference  4  were  used  to  establish  B.  Thus^  it  appears 
that  the  buckling  theory  of  the  present  paper  includirtg  the  modification 
given  in  equation  (8)  will  agree  with  experimental  results  about  as  well 
as  the  methods  available  to  predict  buckling  at  room  temperature  agree 
with  room-temperature  buckling  results.  A  better  evaluation  of  the  theory 
could  be  obtained  from  cylinders  with  large  values  of  r/t  so  that  the 
buckling  temperat\ire  would  be  B.ower  and  in  a  more  usable  temperature  range 
for  the  material. 


CONCLUSIONS 


The  buckling  temperature  of  a  ring-stiffened  cylinder  loaded  in  com¬ 
pression  is  shown  to  be  a  function  of  the  radius -thickness  ratio  of  the 
cylinder^  the  thermal  expansion  coefficient  of  the  cylinder  material,  and 
a  buckling  coefficient  which  is  presented  in  tabular  and  graphical  form 
for  both  clamped  and  simply  supported  cylinders.  The  buckling  coefficient 
is  a  function  of  the  ratio  of  applied  axial  stress  to  the  classical 
buckling  stress  and  the  cylinder  curvature  parameter.  From  the  results 
of  the  analysis,  the  following  conclusions  can  be  draTO: 

1.  For  moderate  to  large  values  of  the  curvature  parameter,  the 
cylinder  buckling  temperature  is  essentially  independent  of  length. 

2.  For  small  values  of  the  radius-thickness  ratio  buckling  of  the 
cylinder  due  to  teraperatuce  alone  occurs  beyond  the  useful  temperature 
range  for  most  materials.  For  larger  values  of  the  radius-thickness  ratio 
such  buckling  can  occur  at  temperatures  within  the  usable  temperature 
range  for  the  material. 


5.  A  cylinder  can  support  a  significant  amount  of  axial  load  without 
reducing  the  buckling  temperature. 

4.  For  moderate  to  large  values  of  the  curvature  parameter,  there 
is  little  difference  between  the  buckling  temperatures  for  clamped  and 
simply  supported  cylinders  of  the  same  proportions  when  the  axial  stress 
approaches  the  classical  buckling  stress.  At  lower  values  of  axial  stress 
there  is  a  substantial  difference. 

5-  The  theoretical  buckling  stress  for  an  unheated  cylinder  under 
uniform  axial  compression  is  about  95  percent  of  the  classical  value  due 
to  the  restraint  of  Poisson  expansion  at  the  ends. 


Ik 


A  modification  of  the  theory  la  given  \rtilch  takes  Into  account  the 
reduction  in  buckling  stress  from  the  classical  value  for  cylinders  under 
axial  compression.  Con^rlson  of  the  modified  theory  vlth  experiment 
indicates  agreement  comparable  to  that  obtained  by  using  existing 
enqpirical  methods  to  predict  room-temperature  coopressive  buckling 
results. 


Langley  Research  Center, 

National  Aeronautics  and  Space  Administration, 
Langley  Air  Force  Base,  Va.,  Jan.  15,  1962. 
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APPENDIX  A 

BUCKLING  OF  A  CIRCULAE  CYLINDER  UNDER  UNIFORM  AXIAL  COMPRESSION 
AND  AN  ARBITRARY  DISTRIBUTION  OF  CIRCUMFERENTIAL  STRESS 


The  equation  of  equilibrium  governing  the  buckling  of  a  cylindrical 
shell  subjected  to  in-plane  direct  stresses  only  is  (see  ref.  5) 


DV^w  +  — 


Et  ,--,-4  ^ 


V 


+  c^t 


+  a^t 

<^x2  Sy2 


=  0 


(Al) 


Dividing  equation  (Al)  by  D  gives 


V^v  +  12  ^  ^  +  kx  —  —  +  k^^  & 


L^  ^  L2  dx2  l2  c)y2  1  + 


n  cQs  ^ 


=  0  (A2) 


where 


Z  =  ilA  - 

rt 


CTX  =  - 


L^t 


and  the  arbitrary  distribution  of  Oy  is  given  by  the  Fourier  series 


P  .  *”  -  c; 

^  - 21-  cos  22^ 

^  L-t  ^  1  +  5on  2 


The  Sn  values  are  defined  so  that  Oy  at  x  equal  zero  is  given  by 


0„(0)  .  ^ 


(A5) 
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This  definition  requires  that 


TO  _ 

H— ^  =  1 

^  1  +  ®0n 


(A4) 


It  can  then  he  seen  that 


is  the  ratio  of  the  amplitude  of  the 


^  ^On 

nth  component  of  stress  to  the  stress  at  x  equal  zero. 


The  method  used  for  the  solution  of  equation  (A2)  is  the  Galerkin 
method.  (See  ref.  J.)  In  this  method  the  differential  equation  of 
equilibrivim  is  represented  as 


h 


Q(v)  =  0 


(A5) 


where  in  equation  (A2) 


ci.v^  +  3:^^r‘‘S!L+i,!ti5L+kjESL7-  _2n_ 


The  deflection  w  is  talen  as 


m=i 

The  set  of  functions  Vjj^  must  satisfy  boundary  conditions  but  not  neces- 
sarilv  the  eauilibrium  equation.  The  coefficients  are  determined 

from  the  condition 


VinQ(w)  dx  dy  =  0 


(A6) 


The  details  of  the  solution  follow 


n  C5\\D  (-■  f 
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Simply  Supported.  Ends 

The  orientation  of  the  coordinate  system  vith  respect  to  the 
cylinder  is  shown  in  figure  12.  An  expression  for  w  which  satisfies 
the  differential  equation  in  the  y-direction  and  the  condition  of  simple 
support  at  the  ends  of  each  bay  is 


w 


sin 


sin 


nmx 

L 


(A?) 


and  Vm  can  be  identified  as 


V„  =  sin  ^  sin  IMX  (AO) 

A  L 


Substitution  of  equations  (A7)  and  (a8)  into  equation  (a6)  and  inte 
grating  over  the  limits  yields 


®ii^V+l 


kp'^ 


XI  ^i  (°in-i 
1=1  ^ 


=m+i) 


=  0 


(m  =  1,  2,  5,  .  .  .)  (A9) 


where  by  definitiOii 


s_i  =  Si 


and 


%  =  |(m  -  1)=  +  ^ 

ii  /  _  V 


n^ljm  -  1)2  +  p2j' 


-  (m  -  l)^kj. 


The  condition  necessary  for  buckling  is  that  the  determinant  of  the 
coefficients  of  equations  (A9)  eqvial  zero.  This  determinant  can  be 
written  as 
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» 


n 


1 


“1 


£Mg 


®0  ®2 


“2 

+  Sj 


“5 

-Bg  +  Sl| 


m  =  2 


-"l  +  «3 


2M:. 

-80+64  -Si  + 


m  =  5 


-Bg  +  84 


a«4 

-®T  +  ^  -  “0  +  =6 

p'^k 


* 


=■  0  (AlO) 


For  kjj  eq.ual  zero,  this  result  corresponds  to  that  obtained  in  refer¬ 
ence  2  for  buckling  of  a  simply  supported  cylinder  subjected  to  an 
arbitrary  distribution  of  circumferential  stress.  Reference  2  is  based 
on  Donnell's  eighth-order  differential  equation. 

Equation  (AlO)  may  be  solved  for  k  for  given  values  of  k^  if 

the  Sfj  values  are  knovn.  By  using  an  increasing  number  of  teiros  and 
equations,  the  solution  vill  converge.  The  solution  for  k  should  be 
minimized  vith  respect  to  (i  to  give  the  lowest  buckling  stress.  The 
values  of  p  should  be  restricted  so  that  an  integral  number  of  waves 
appear  around  the  cylinder  circumference.  However,  this  number  is 
usually  large  enough  that  p  may  be  thought  of  as  a  continuous  function. 


is 


Clamped  Ends 

An  expression  for  w  which  satisfies  the  condition  of  clamped  ends 


=  sin  ^  % 

^  m=l 


cos(m  -  1)—  -  cos(m  +  l)  — 
L  L 


(All) 


and  Vjjj  is  therefore 


''m  "  =1”  f 


cos(m  -  1)—  _  cos(m  +  1)— 
L  L 


(A12) 
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Performing  a  similar  operation  to  that  for  the  simply  supported  cylinder 
yields 


^Im)  ~  ®’m-2Mni(l  ~  ^Im)  (l  ~  ^2m)  “  ^+2^+2* 

2  00 

H  (sm-i-2  -  2sin_i  +  “  SnH-i-2  +  = 

1=1  '  ' 


(m  =  1,  2,  3,  .  .  .)  (A13) 


where 


Setting  the  determinant  of  the  coefficients  of  equations  (AI5)  equal  to 
zero  yields  the  buckling  stress 


''1 


1  p(' 


4-  K; 


I 

:  j  +  Jnj  -  r.^ 


-  Jcq  +  -?;>  -  c;*  -2r>j  ♦  Jr.j  *  s., 


+■  Cq  -  352  Jsi*  -• 


m  »  3 


•Sin 

♦  so  - 


ySp  4-  5^1*  *  t'.. 


•  Cca  +  Tn  +  Csi.  •  }u: 


-s^  +  Cc^  -  r-. 


P"-k 


P*-)^ 


i  -  0  (AlU) 


2  -  s;^  +  20^  .  SQ 


The  solution  of  equation  (A14)  is  similar  to  that  indicated  for  equa¬ 
tion  (AIO)  . 
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APPENDIX  B 

STRUCTURAL  BEHAVIOR  OF  RING-STIFFENED  CYLINDER  LOADED  IN 
AXIAL  COMPRESSION  AND  HEATED  UNIFORMLY 


Thermal  Stress 

It  is  ass\imed  that  the  skin  is  at  a  constant  temperature  and  the 
rings  are  at  some  lower  constant  temperature^  the  difference  in  tempera¬ 
ture  heing  AT.  .Also  there  are  many  rings  so  that  each  hay  may  he 
assumed  to  have  the  same  deflection  curve  which  will  he  symmetrical 
about  the  middle  of  the  hay.  Deflections  will  not  vary  around  the 
circumference  and  will  he  given  hy  the  solution  of  the  following  (see 
ref .  6) : 


D  ^  +  a.t  ^  +  -5^-  Et 
dx^  ^  dx2  r2 


^  x=0  -  V  -  (Bl) 

x=L  [ 


dw 

^  x=0 
x=L 


=  0 


Deflections  are  measui'ed  positive  inward  from  the  position  an  -unstressed 
cylinder  would  assume  without  rings  and  subjected  to  a  uniform  tempera¬ 
ture  T.  The  solution  for  w  is 


.00  M 

0  sinh  —  cosh  — (-9  sin  ^  cos 
2  2  2 


—r  10  sinh  —  cos  ^ 

0  [y  2  2 
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where 


0  =  Z  |/l  -  7 

0  =Jl|5  Z  /l  +  7 


^/5  z 


kqi  — 


-  <^’^n 

rt2 


12(1 


The  deflection  v  will  cause  a  circumferential  strain  pj  thus^  hy 
will  he  given  from  the  stress-strain  relations  as 


ff 


y 


=  P<Jx  + 


yfs 

r 


(B5) 


In  order  to  apply  this  result  to  the  huckling  analysis,  It  Is  necessary 
to  expand  a-y  In  a  J'OTxrier  series.  This  expansion  can  he  done  hy 
substituting  equation  (B2)  into  equation  (B5)  and  expanding  this  result 
In  a  series  recalling  the  defining  relationship  given  in  equations  (AJ) 
and  (a4-)  .  The  resulting  Poirrier  coefficients  s^j  are  given  hy 


P^adf  M\(q2  j.  m2) 


^02  02  ^  _  402jj2j{2  ^0  sinh  0+9  sin  0^ 


Sn  -  0 


even/ 


(n,  odd) 


>  (34) 


and  k  Is  identified  as 


(B5) 


k  =  (ikjj  +  kj 
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The  only  stress  component  which  produces  a  net  force  over  the  cylindei’ 


length  is  sq;  thus 


^0 

—  equals  — the  ratio  of  the  average  cir- 
2  ^yL^' 


cumfermtial  stress  to  the  stress  at  a  ring. 


Buckling 


After  the  thermal  stress  is  determined^  the  results  can  be  used  in 
the  buckling  analysis  of  appendix  A  to  calculate  the  buckling  tempera¬ 
ture.  Inasmuch  as  the  odd  numbered  s^  are  zero,  each  stability  deter¬ 
minant  (eqs.  (AIO)  and  (AlU) )  can  be  factored  into  two  determinants,  one  C 

containing  terms  involving  the  odd-numbered  deflection  coefficients  and  * 

txie  other  containing  terms  involving  even-numbered  deflection  coeffi¬ 
cients  .  The  two  determinants  correspond  to  deflection  patterns  symmet¬ 
rical  and  anti symmetrical  about  a  plane  perpendicular  to  the  axis  of  the  » 

cylinder  at  the  midlength  point.  The  two  determinants  for  the  simply 
supported  cylinders  are 


m  =  1 


a  =  5 


ad> 


-  Sg  +  Sj 


^3 

+  sjj 


-“It  +  =6 


JN., 


-=it  +  “6 


P*"  i  ***  ^ 

-=2  +  “3 


‘0  -''6 


-E2  +  S3 


=  0  (b6) 


+  loj 


Co  + 


m  =  2 


2M, 


I  p2(ttk,  + 


-  -  to  “i; 


RU 

-S2  > 


“2  +  =6 


-=it  +  “3 


aitj 


+  1«T,' 

-“2  +  “10 


_  -  so  +  SQ 


“6 

-E4  +  E3 

-So  +  s^o 


''T, 


■'0  ■'12  •  • 


=  0  (B7) 
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For  the  clamped  cylinder  the  following  two  determinants  are 
obtained: 


Values  of  calculated  by  using  the  various  stability  determi¬ 

nants  revealed  a  more  suitable  buckling  parameter.  The  new  i)arameter  t 
is  such  that  the  buckling  temperature  is  given  in  terms  of  r/t  rather 
than  L/t  as  indicated  in  the  following  equation: 


+ 


T  t 
a  r 


(BIO) 


The  parameter  t  is  similar  to  the  parameter  7  as  kp  is  similar  to 
k^.  The  parameter  t  is  related  to  kp  by 


(Bll) 


Values  of  kp  were  found  to  range  through  several  orders  of  magnitude 

with  increasing  values  of  Z  whereas  x  shows  only  a  small  variation 
as  shown  in  figures  1  and  2.  Hence,  results  are  better  shown  in  terms 
of  T  and  values  are  more  easily  interpolated. 
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Nonrigid  rings 

In  order  to  correct  for  ring  flexibility  it  vill  be  assiimed  that 
ring  deflections  can  be  neglected  dxiring  buckling^  but  the  effect  of 
these  defle  -uions  will  be  included  in  calculating  the  circumferential 
thermal  stresses.  In  other  words^  the  ring  is  assumed  to  deflect' 
because  of  the  application  of  and  by  stresses  but  is  assumed  to 

be  rigid  during  the  buckling  process.  This  correction  is  accomplished 
as  follows.  The  first  boundary  condition  on  w  in  equation  (B1) 
becomes 


w 


x=0 

x=L 


(B12) 


where  '4f  is  the  deflectional  stiffness  of  the  ring  and  by  is  the 
average  circumferential  stress  in  the  skin.  For  a  ring  that  resists 
radial  deflection  by  stretching  in  the  circumferential  direction,  t 
is  given  by 


(B15) 


With  the  boundary  condition  (B12)  and  repeating  the  thermal  stress 
analysis  as  before,  it  is  found  that  k  can  be  taken  as 


.  T 

k  =  pkjj  +  kjp 


12(1  -  -  ELt 


(B14) 


with  Sj^  defined  as  before  and  k<r 
cient  for  nonrigid  rings.  However, 


the  buckling  temperature  coeffi- 
by  is  related  to  and  sq  by 


^0  (t\^ 

^  12(1  - 


(BI5) 


and  k^  is  therefore  given  by 


,  -  ,  Sq  ELt\  , 


2  j.2^ 


(BI6) 


1-1  r-|  ON  VO  O 
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This  equation  can  he  rewritten  in  terms  of  t',  t,  and  y,  and  with 
the  ring  stiffness  given  "by  equation  (BIJ),  equation  (Bl6)  appears 
as 


=  T  + 


=■  +  T 


5(1  -  u2) 


(B17) 


Thus,  the  buckling  temperature  coefficient  for  cylinders  with  flexible 
rings  equals  the  coefficient  for  rigid  rings  plus  an  additional  term 
which  is  proportional  to  the  average  circumferential  stress  in  the 
cylinder . 
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Figure  1.-  Temperature  buckllne  coefficient  for  axially  compressed 
simply  supported  ring-stiffened  cylinders. 


Figure  2.-  Temperature  buckling  coefficient  for  axially  compressed 
clamped  ring- stiffened  cylinders. 
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Figure  5.-  Plot  Illustrating  change  of  buckle  shape  as  axial  stress 

is  Increased,  p  =  h/'h. 
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Figure  7*~  Comparison  of  the  average  thermal  "buckling  stress  with  the 
"buckling  stress  for  rmiform  circumferential  compression. 
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Figure  8.-  Ratio  of  average  clrcxnnferential  stress  to  stress  at  a  ring 
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